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Abstract 

A possibility to measure sin^ 2^13 using reactor neutrinos is examined in detail. It is shown that 
the sensitivity sin^ 2^13 > 0.02 can be reached with 40 ton-year data by placing identical CHOOZ- 
like detectors at near and far distances from a giant nuclear power plant whose total thermal energy 
is 24.3 GWth- It is emphasized that this measurement is free from the parameter degeneracies which 
occur in accelerator appearance experiments, and therefore the reactor measurement plays a role 
complementary to accelerator experiments. It is also shown that the reactor measurement may be 
able to resolve the degeneracy in 023 if sin^ 20i3 and cos^ 2023 ^-^e relatively large. 
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I. INTRODUCTION 


Despite the accumulating knowleciges of neutrino masses and the lepton flavor mixing 
by the atmospheric the solar j^, 1 ^, and the accelerator neutrino experiments, the 
(1-3) sector of the Maki-Nakagawa-Sakata (MNS) matrix is still in the dark. At the 
moment, we only know that \Ue^\ = sin^is = S 13 is small, < 0.03, by the bound imposed 


by the CHOOZ reactor experiment p. In this paper we assume that the light neutrino 
sector consists of three active neutrinos only. One of the challenging goals in an attempt to 
explore the full structure of lepton flavor mixing would be measuring the leptonic CP or T 
violating phase 5 in the MNS matrix. If KamLAND confirms the Large-Mixing-Angle 
(LMA) Mikheev-Smirnov-Wolfenstein (MSW) [1, solution of the solar neutrino problem, 
the most favored one by the recent analyses of solar neutrino data BII 3 , we will have 
an open route toward the goal. Yet, there might still exist the last impasse, namely the 
possibility of too small value of ^ig. Thus, it is recently emphasized more and more strongly 
that the crucial next step toward the goal would be the determination of ^ig. 

In this paper, we raise the possibility that Pg disappearance experiment using reactor 
neutrinos could be potentially the fastest (and the cheapest) way to detect the effects of 
nonzero 6 ^ 13 . In fact, such an experiment using the Krasnoyarsk reactor complex has been 
described earlier in which the sensitivity to sin^ 2^13 can be as low as ~ 0 . 01 , an order 
of magnitude lower than the CHOOZ experiment. We will also briefly outline basic features 
of our proposal, and reexamine the sensitivity to sin^ 2^13 in this paper. 

It appears that the most popular way of measuri^^ig is the next generation long baseline 
(LBL) neutrino oscillation experiments, MINOS [l^, OPERA 3. and the JHF phase I 
l3. It may be followed either by conventional superbeam m experiments, the JHF phase 
■ I |l4| and possibly others ISIS, or by neutrino factories [l^ la. t is pointed out, 
however, that the measurement of ^13 in LBL experiments with only neutrino channel (as 
planned in the JHF phase I) would suffer from large intrinsic uncertainties, on top of the 
experimental errors, due to the dependence on an unknown CP phase and the sign of 
|2flj |. Furthermore, it is noticed that the ambiguity remains in determination of ^13 and 
other parameters even if precise measurements of appearance probabilities in neutrino as 
well as antineutrino channels are carried out, the problem of the parameter degeneracy 
0 , H H HI H 0 0 . (For a global overview of the parameter degeneracy, see 0 .) 
While some ideas toward a solution are proposed the problem is hard to solve experimentally 
and it is not likely to be resolved in the near future. 

We emphasize in this paper that reactor Ve disappearance experiment provide particularly 
clean environment for the measurement of ^ 13 . Namely, it can be regarded as a dedicated 
experiment for determination of 6 ^ 13 ; it is insensitive to the ambiguity due to all the remaining 
oscillation parameters as well as to the matter effect. This is in sharp contrast with the 
features of LBL experiments described above. Thus, the reactor measurement of ^13 will 
provide us valuable information complementary to the one from LBL experiments and will 
play an important role in resolving the problem of the parameter degeneracy. It will be 
shown that reducing the systematic errors is crucial for the reactor measurement of 6^13 to 
be competitive in accuracy with LBL experiments. We will present a preliminary analysis 
of its possible roles in this context. It is then natural to think about the possibility that one 
has better control by combining the two complementary way of measuring 6 * 13 , the reactor 
and the accelerator methods. In fact, we will show in this paper that nontrivial relations 
exist between the ^13 measurements by both methods thanks to the complementary nature 
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of these two methods, so that in the luckiest case one may be able to derive constraints on 
the value of the CP violating phase 6, or to determine the neutrino mass hierarchy. 


II. REACTOR EXPERIMENT AS A CLEAN LABORATORY FOR 613 MEA¬ 
SUREMENT 


Let us examine in this section how clean the measurement of 6*13 by the reactor experi¬ 


ments is. To dehne our notations, we note that the standard notation j27 


U = 


C12C13 S12C13 Si3e“ 

—-S12C23 — Ci 2S23S13C‘^ C12C23 — Sl2S23Sl3e*'^ S23C13 , 

512-523 ~ Ci2C23Si3e*'^ —Ci2523 “ 5i2C23Si3e®'^ C23C13 


( 1 ) 


is used for the MNS matrix throughout this paper where Cij and {i,j = 1 — 3) imply 
cos 9ij and sin 9ij, respectively. The mass squared difference of neutrinos is dehned as Am^^- = 
where is the mass of the zth eigenstate. 

We examine possible ’’contamination” by 6, the matter effect, the sign of Am|^, and the 
solar parameters one by one. We hrst note that, due to its low neutrino energy of a few 
MeV, the reactor experiments are inherently disappearance experiments, which can measure 
only the survival probability P(he ^e)- It is well known that the survival probability does 
not depend on the CP phase 6 in arbitrary matter densities | 2 ^ . 

In any reactor experiment on the Earth, short or long baseline, the matter effect is very 
small because the energy is quite low and can be ignored to a good approximation. It can 
be seen by comparing the matter and the vacuum effects (as the matter correction comes in 
only through this combination in the approximate formula in 
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where 


Ajj — 


Am?L 




2E 


( 3 ) 


with E being the neutrino energy and L baseline length. The best £t value of |Am|]^| is 
give n by |Am|]^| = 2.5 x lO-^eV^ from the Super-Kamiokande atmospheric neutrino data 
j 20 | . and we use this as the reference value for jAmg^l throughout this paper, a = v^G^iVe 
denotes the index of refraction in matter with Gp being the Fermi constant and W the 
electron number density in the Earth which is related to the Earth matter density p as 
Nf. = Yepjrrtp where W is proton fraction. Once we know that the matter effect is negligible 
we immediately recognize that the survival probability is independent of the sign of Amg^. 

Therefore, the vacuum probability formula applies. The general probability formula in 
vacuum is analytically written as j27| 


f P(n„ ^ Up) 1 

1 ^p) f 



( 4 ) 
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where a,f 3 = e,/i, r, and the minus and signs in front of the Im (UajU^jU*j.Ui3k) term in 
0 correspond to neutrino and antineutrino channels, respectively. From © the exact 
expression for P{ 9 e ^e) is given by 

1 - P{Ue ^ Pe) = 4 ^ \Uej 
j<k 

= sin^ 26*13 sin^ 26*13 sin A31 sin A21 

+ (ci3 sin^ 26*12 + Ciasin^ 26*13 cos A31) sin^ ( 5 ) 

where the parametrization m has been used in the second line. The three terms in the 
second line of © are suppressed relative to the main depletion term, the hrst term of 
the right-hand-side of © , by e, e^/sin^ 26*13, eT respectively, where e = Am2i/|Am|;^|. 
Assuming that lAm^J = ( 1 . 6 - 3 . 9 ) x 10 “^eV^ ji^, e ~ 0 . 1 - 0.01 for the LMA MSW solar 
neutrino solution |^, Tlol | . Then, the hrst and the third terms in the second line can be 
ignored, although the second term can be of order unity compared with the main depletion 
term provided that e ~ 0 . 1 . (Notice that we are considering the measurement of sin^ 26*13 in 
the range of 0 . 1 - 0 . 01 .) Therefore, assuming that |Am|i| is determined by LBL experiments 
with good accuracy, the reactor Pg disappearance experiment gives us a clean measurement 
of 6*13 which is independent of any solar parameters except for the case of high Am^i LMA 
solutions. 

If the high Am^i LMA solution with Am^i ~ 10 “^ eV^ turns out to be the right one, 
we need a special care for the second term of the second line of ©. In this case, the 
determination of ^13 and the solar angle ^12 is inherently coupled,^ and we would need joint 
analysis of near-far detector complex (see the next section) and KamLAND. 
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III. NEAR-FAR DETECTOR COMPLEX: BASIC CONCEPTS AND ESTIMA¬ 
TION OF SENSITIVITY 

In order to obtain good sensitivity to sin^ 26*13, selection of an optimized baseline and 
having the small statistical and systematic errors are crucial. For instance, the baseline 
length that gives the oscillation maximum for reactor Pg’s which have typical energy 4 MeV 
is 1 . 7 km for Am^ ~ 2.5 x 10 “^ eV^. Along with this baseline selection, if systematic and 
statistical errors can be reduced to 1 % level, which is 2.8 times better than the CHOOZ 
experiment , an order of magnitude improvement for the sin^ 2^13 sensitivity is possible 
at Am^ ~ 2.5 x 10 “^ eV^. In this section we demonstrate that such kind of experiment 
is potentially possible if we place a CHOOZ-like detector at a baseline 1.7 km in 200 m 
underground near a reactor of 24.3 GWth thermal power. The reactor can be regarded as 
a simplihed one of the Kashiwazaki-Kariwa nuclear power plant which consists of seven 
reactors and whose maximum energy generation is 24.3 GWth- 

Major part of systematic errors is caused by uncertainties of the neutrino flux calculation, 
number of protons, and the detection efficiency. For instance, in the GHOOZ experiment, 
the uncertainty of the neutrino flux is 2.1 %, that of number of protons is 0.8 %, and that of 


^ The effect of nonzero 6*13 for measurement of 612 at KamLAND is discussed in 
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detection efficiency 1.5% as is shown in the Table IH The uncertainty of the neutrino flux 
includes ambiguities of the reactor thermal power generation, the reactor fuel component, 
the neutrino spectra from hssions, and so on. The uncertainty of the detection efficiency 
includes systematic shift of dehning the hducial volume. These systematic uncertainties, 
however, cancel out if identical detectors are placed near and far from the reactors and data 
taken at the detectors are compared.^ 

To estimate how good the cancellation will be, we study the case of the Bugey experiment, 
which uses three identical detectors to detect reactor neutrinos at 14/40/90 m. For the Bugey 
case, the uncertainty of the neutrino flux improved from 3.5 % to 1.7% and the error on the 
solid angle remained the same (0.5% —*• 0.5%). If each ratio of the improvement for the 
Bugey case is directly applicable to our case, the systematic uncertainty will improve from 
2.7% to 0.8% as shown in the Table HI The ambiguity of the solid angle will be negligibly 
small because the absolute baseline is much longer than the Bugey case. We are thinking 
of a case that a front detector is located at 300 m away from the reactor we consider. In 
the actual setting with the Kashiwazaki-Kariwa power plant two near detectors may be 
necessary due to extended array of seven reactors. Hereafter, we take 2% and 0.8% as the 
reference values of the relative systematic error agys for the total number of z/g events in 
our analysis. Let us examine the physics potential of such a reactor experiment assuming 
these reference values for the systematic error. We take, for concreteness, the Kashiwazaki- 
Kariwa reactor of 24.3 GWth thermal power and assume its operation with 80 % efficiency. 
Two identical liquid scintillation detectors are located at 300 m and 1.7 km away from the 
reactor and assumed to detect Ue by delayed coincidence with 70 % detection efficiency. The 
UeS of 1-8 MeV visible energy, — 0.8 MeV, are used and the number of events are 

counted in 14 bins of 0.5 MeV. Without oscillation, a 10 (40) ton-year measurement at the 
far detector yields 20,000 (80,000) hg events which is naively comparable to a 0.7 (0.35) % 
statistical error. 

First, let us calculate how much we could constrain sin^ 20 i 3 . Unlike the analysis in 
which uses the ratio of the numbers of events at the near and the far detectors, we use the 
difference of the numbers of events Ni(Lo) — [Li/L 2 )‘^Ni{Li), because the statistical analysis 
with ratios is complicated. (See, e.g., jS^-) The dehnition of Ay^, which stands for the 
deviation from the best £t point (non-oscillation point) is given by 


Ay2(sin2 2ffi3, \Aml^\) 

14 I ^i(0)(-b2) - Ai(o)(Ti) 

i=l 


v(i2) - (fe) A.Ui 


V(„)(i2) + (It) V(„)(il) + 


( 6 ) 


Ni{Lj) = W(sin^26'i3, \Aml^\]Lj), Ai(o)(Tj) = Ai(0,0;Lj), 

where is the relative systematic error for each bin which is assumed to be the same 
for all bins and Aj(sin^ 2 ^ 13 , lAml^l) denotes the theoretical number of Pg events within the 


^ This is more or less the strategy taken in the Bugey experiment m- The Krasnoyarsk group also plans 
in their Kr2Det proposal to construct two identical 50 ton liquid scintillators at 1100 m and 150 m 
away from the Krasnoyarsk reactor. They indicate that the systematic error can be reduced to 0.5% by 
comparing the front and far detector. 
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ith energy bin. In principle both the systematic errors cr^b” sys (absolute normalization) and 
(relative normalization) appear in the denominator of (jHI), but by taking the difference, 
we have (1 + a^(^.,yj[(l + <^)iV,(L 2 ) - {U/L,fN,{L,)] - [iV^L^) - {L,/L,fN,{U)] = 
(j^ygiVi(L 2 ) + crab°sys[^*(-^ 2 ) “ {Li/L 2 YNi{Li)\ which indicates that the systematic error 
is dominated by the relative error as the second term [Ni{L 2 ) — (Li/L 2 )^iVj(Li)] is 
supposed to be small. In fact we have explicitly verihed numerically that the presence of 
(^abssys)^[-^*(-^ 2 ) “ {Li/L 2 YNi{Liyf‘ in the denominator of (0 does not affect any result. 
From the assumption that the relative systematic error for each bin is distributed equally 
into bins, is estimated from the relative systematic error crgys for the total number of 
events by 
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(y“(i2))= 


y“(i2) = i]iv.,„,(L2), 

i 


(7) 


since the uncertainty squared of the total number of events is obtained by adding up the bin- 
by-bin systematic errors ((T^yg)^iVj(Q)(L 2 ); The ratio cTgy“/(Tsys is about 3 in our analysis. In 
Fig. in the 90 % CL exclusion limits, which corresponds to = 2.7 for 1 degree of freedom, 
are presented for two cases: a 10 ton-year measurement with the 2 % systematic error of the 
total number of events and a 40 ton-year measurement with the 0.8% error. The hgure 
shows that it is possible to measure sin^ 26^13 down to 0.02 at the maximum sensitivity with 
respect to |Am|]^|, and to 0.04 for larger |Am|^| by a 40 ton-year measurement, provided the 
quoted values of the systematic errors are realized. The CHOOZ result is also depicted 
in Fig. n For a fair comparison with the CHOOZ contour, we also present in Fig. Q the 
results of analysis with 2 degrees of freedom, which correspond to Ay;^ = 4.6 for 90% CL, 
without assuming any precise knowledges on |Am|]^|. 

Next, let us examine how precisely we could measure sin^26'i3. The dehnition of Ay^ is 


Ax^(sin2 20i3, lAm^il) 

14 I iV,(best)(L2)- (t^)'iV,(best)(Tl) - A, (L2) - ' A, (Li) 

^i(best)(F2) + Aj(best)(Tl) + (o'^ys (-^2) 


( 8 ) 


where Aj(best) denotes Aj for the set of the best £t parameters (sin^ given 

artihcially. ( 73 “ is obtained in dZD by replacing Aj(o) with Ai(best) and the ratio cr^yg/cTsys is 
about 3 again. We assume that the value of jAml^l is known to a precision of 10“^ eV^ 
by the JHF phase I by the time the reactor measurement is actually utilized to solve the 
degeneracy. Then, we rely on the analysis with 1 degree of freedom, hxing lAml^l as 
= 2.5 X 10“^ eV^. The 90% CL allowed regions of 1 degree of freedom, whose 
bounds correspond to Ay^ = 2.7, are presented in Fig. |21 for the values of sin^ from 

0.05 to 0.08 (0.02 to 0.08) in the unit of 0.01 in the case of a 10 ton-year (40 ton-year) 
measurement with systematic error cisys = 2.0(0. 8 ) %. We can read off the errors at 90 % CL 
in sin^ 26^13 and it is almost independent of the central value sin^ 2013*^^*^ Thus, we have 

sin^26'i3 = sin^ 26*[3^^*^ ± 0.043 (at 90% CL, d.o.f. = 1) 
for sin^ 26'[3®^*^ > 0.05 
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in the case of dgys = 2 % with a 10 ton-year measurement, and 

sin^ 26*13 = sin^ 26 * 13 ®'^*^ ± 0.018 (at 90 % CL, d.o.f. = 1 ) 
for sin^ > 0.02 

in the case of agys = 0.8 % with a 40 ton-year measurement. 


IV. THE PROBLEM OF THE (6*13, 6 , 6*23, Am|i) PARAMETER DEGENERACY 


We explore in this and the following sections the possible signihcance of reactor measure¬ 
ments of 6*13 in the context of the problem of the parameter degeneracy. We show that the 
reactor measurement of 6*13 can resolve the degeneracy at least partly if the measurement is 
sufficiently accurate. Toward the goal we hrst explain what is the problem of the parameter 
degeneracy in long-baseline neutrino oscillation experiments. It is a notorious problem; a set 
of measurements of the disappearance probability and the appearance oscillation proba¬ 
bilities of z/g and —>• Ug, no matter how accurate they may be, does not allow unique 
determination of 6*13, S, and 6*23. The problem was hrst recognized in the form of intrinsic 
degeneracy between the two sets of solutions of (6*93 , 6*13) for a given set of measurements 


in two different channels z/,. 


zz, 


e and 


It was then observed independently 


that the similar degeneracy of solutions of (6*13, 5 ) exists in measurement of Ug appearance 
in neutrino and antineutrino channels | 2 ^. They made the hrst systematic analysis of the 
degeneracy problem. It was noticed that the degeneracy is further duplicated provided that 
the two neutrino ma ss p atterns, the normal (Am^^ > 0) and the inverted (Am^^ < 0) hier¬ 
archies, are allowed j 23 |. Finally, it was pointed out that the degeneracy can be maximally 


setting in 


26 | 


eight-fold m- Analytic structure of the degenerate solutions was worked out in a general 


To illuminate the point, let us hrst restrict our treatment to a relatively short baseline 
experiment such as the CERN-Frejus project In this case, one can use the vacuum 

oscillation approximation for the disappearance and the appearance probabilities. From the 
general formula 0) we have 


1 - P{u^ = 


4 AA-t 

j<k 

A31 

23 ■ ‘ 


AE 


= sin^ 26*23 sin^ 


—C32 sin^ 26*23 — 513523 sin 26*23 sin 26*12 cos 6 ) sin A21 sin A31 

± 0 ( 6 ^) ± 0 ( 4 ), 


( 9 ) 


7 









u^) 


( A™2 t 

' /\m? L 

T2^Im(f/„j(7;^(7if/.t)sm' 

j<k 

2 • 2 ofl -2 ^31 

S 23 Sin 26^13 Sin —^ 


2 E 


1 


Jr sin A21 sin A31 cos 6 


A3 


=F Jr sin A21 sin"^ sin (5 + 0(6833), 


( 10 ) 


where e = Am^^/IAml^l, Jr = sin 26*23 sin 2012C33S13, cind the parametrization has been 
used in the second line in each formula. The minus and plus signs in front of sin 5 term in 
(uni) correspond to neutrino and antineutrino channels, respectively. An explicit perturbative 
computation in indicates that the matter effect enters into the expression in a particular 
combination with other quatities (in the form of s\^aL /so that the effect is small. By 
the disappearance measurement at JHF, for example, sin^ 26* 93 a nd jAml^l will be determined 
with accuracies of 1 % for 0.92 < sin^ 26*23 A 1-0 (Fig. 11 in [l^), and 4 % level, respectively 
l 3 |.^ If 6*23 is not maximal, then we have two solutions for 6*23 (^23 and n /2 — 6*23), even if 

we ignore the uncertainty in the determination of sin^ 26*23. For example, if sin^ 26*23 = 0 . 95 , 

2 

23 


which is perfectly allowed by the most recent atmospheric neutrino data | 29 | 1 , then s;, can 


be either 0.39 or 0 . 61 . Since the dominant term in the appearance probability depends upon 
S23 instead of sin^2023, it leads to ±20% difference in the number of appearance events in 
this case. On the other hand, in the case of maximal mixing, it still leaves a rather wide 
range of 6*23, despite such fantastic accuracy of the measurement. 1 % accuracy in sin^ 26*23 
implies about 10 % uncertainty in 833. Thus, whenever we try to determine sin^ 26*13 from 
the appearance measurement, we have to face the ambiguity due to the two-fold nature of 
the solution for 833. 

Let us discuss the simplest possible case, the LOW or the vacuum (VAC) oscillation 
solution of the solar neutrino problem. (See e.g., for a recent discussion.) In this case, one 
can safely ignore terms of order e in Q and m- Then we are left with only the hrst terms 
in the right-hand-side of these equations, the one-mass scale dominant vacuum oscillation 
probabilities. Now let us dehne the symbols x = sin^ 26*13 and y = 833. Then, Q and (ITn|l 
take the forms y = yi or j/2 (corresponding to two solutions of 833) and xy = constant, 
respectively, for given values of the probabilities. It is then obvious that there are two 
crossing points of these curves. This is the simplest version of the (6*13,^23) degeneracy 
problem. We next discuss what happens if e is not negligible though small: the case of LMA 
solar neutrino solution. In this case, the appearance curve, xy = constant, split into two 
curves (though they are in fact connected at their maximum value of 8^3) because of the 
two degenerate solution of the set ( 5 , 6*13) that is allowed for a given set of values of 8^3, 
P{^ii ^e) and P{nfj, —>■ he)- Then, we have, in general, four crossing points on the x-y plane 
for a given value of sin^ 26*23, the four-fold degeneracy. Simultaneously, the two y = constant 


^ Usually one thinks of determining not but |Am| 2 | by the disappearance measurement. But, it 

does not appear possible to resolve difference between these two quantities because one has to achieve 
resolution of order e for the reconstructed neutrino energy. 
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lines are slightly tilted and the splitting between two curves becomes larger at larger sin^ 26 * 13 , 
though the effect is too tiny to be clearly seen. If the baseline distance is longer, the Earth 
matter effect comes in and further splits each appearance contour into two, depending upon 
the sign of Am|^. Then, we have four curves (or, two continuous contours each of which 
intersects with y = constant line twice) and hence there are eight solutions as displayed 
in Fig. 121“^ This is a simple pictorial representation of the maximal eight-fold parameter 


degeneracy [ 2 ^. To draw Fig.|21 we have calculated disappearance and appearance contours 
by using the approximate formula derived by Cervera et al. H. We take the baseline 
distance and neutrino energy as L = 295 km and E = 400 MeV with possible relevance to 
JHF project 3|. The Earth matter density is taken to be p = 2.3 g ■ cm“^ based on 
the estimate given in The electron fraction W is taken to be 0.5. We assume, for 

dehniteness, that a long-baseline disappearance measurement has resulted in sin^ 26*23 = 0.92 

For the LMA solar neutrino parameters we take tan^ 6*12 = 0.38 
js^. We take the values of these parameters and the matter 
density throughout this paper unless otherwise stated. The qualitative features of the hgure 
remain unchanged even if we employ values of the parameters obtained by other analyses. 


and = 2.5 x 10 ^eV^. 
and Am^i = 6.9 x 10“®eV^ 


V. RESOLVING THE PARAMETER DEGENERACY BY REACTOR MEA¬ 
SUREMENT OF 6*13 


Now we discuss how reactor experiments can contribute to resolve the parameter degen¬ 
eracy. To make our discussion as concrete as possible we use the particular long-baseline 

to illuminate the complementary role played by re- 
It is likely that the experiment will be carried out 
for a number of reasons: the dip in 


experiment, the JHF experiment [14 


actor and long-baseline experiments, 
at around the hrst oscillation maximum (IA 31 I = vr 
energy spectrum in disappearance channel is the deepest, the number of appearance events 


are nearly maximal 
each mass hierarchy 



and the two-fold degeneracy in 5 becomes simple (J <->■ vr — J) for 
1^.^ With the distance L = 295km, the oscillation maximum is at 


around E = 600 MeV. We take the same mixing parameters as those used in Fig. [Hj 


A. Illustration of how reactor measurement helps resolve the (^13, ^23) degeneracy 

Let us hrst give an illustrative example showing how reactor experiments could help 
resolve the (6*13, 6*23) degeneracy. To present a clear step-by-step explanation of the relation¬ 
ship between LBL and reactor experiments, we hrst plot in Fig. |3]the allowed regions in the 
sin^ 26 *i 3 -s |3 plane by measurements of —>■ r'e) alone and —>■ Ve) alone separately. 
The former is indicated by the regions bounded by black lines and the latter by gray lines. 
The solid and dashed lines are used for cases with positive and negative Am^^. The values 
of disappearance and appearance probabilities are chosen arbitrarily for illustrative purpose 


^ The readers might be curious about the feature that the two contours are connected with each other at 
a large 5^3 point. Because J is a phase variable, the contours must be closed as 5 varies. 

® In order to have this reduction, one has to actually tune the energy spectrum so that cos 5 term in m 
averaged over the energy with the neutrino flux times the cross section vanishes, which is shown to be 
doable in Q 
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and are given in the caption of Fig. 0 ] Notice that the negative cnrve is located right 

(left) to the positive curve in neutrino (antineutrino) channel. A plot with only mea¬ 

surement in neutrino mode goes beyond academic interest because the JHF experiment is 
expected to run only with the neutrino mode in its hrst phase. We observe that there is large 
intrinsic uncertainty in the ^13 determination due to unknown 6, the problem addressed in 
[2^. The two regions corresponding to positive and negative Am|]^ heavily overlap due to 
small matter effect. When two measurements of u and u channels are combined, the allowed 
solution becomes a line which lies inside of the overlap of the u and z/ regions for each sign of 
Am|^ in Fig.| 3 ] ® In Fig.^lwe have plotted such solutions as two lines, one for positive Am|^ 
(the solid curve) and the other for negative Amgj^ (the dashed curve) at the hrst oscillation 
maximum IA31I = tt. It may appear curious that the two curves with positive and negative 
Arug^ almost overlap with each other in Fig. El In fact, a slight splitting between the solid 
(Am|g > 0) and dashed (Am|g < 0) lines is due to the fact that both e and the matter 
effect in the case of the JHF experiment are small. Thus, the degeneracy in the set (6*13, 
6*23) is effectively two-fold in this case. 

To have a feeling on whether the reactor experiment described in Sec. Ill will be able 
to resolve the degeneracy, we plot in Fig. El two sets of degenerate solutions by taking a 
particular value of 6*23, sin^ 26*23 = 0 . 92 , the lower end of the region allowed by Super- 
Kamiokande. We denote the true and fake solutions as (sin^ 26*13, S23) and (sin^ 26*^3, s^g^), 
respectively, assuming the true 6*23 satishes 6*23 < 7r/4. We overlay in Fig. El a shadowed 
region to indicate the accuracy to be achieved by the reactor measurement of ^13. If the 
experimental error ( 5 re(sin^ 26*13) ia the reactor measurement of sin^20i3 is smaller than the 
difference 


( 5 de(sin^ 26*13) = I sin^ 26*^3 — sin^ 26*13] (11) 

due to the (6*13, 623) degeneracy, then the reactor experiment may resolve the degeneracy. 
Notice that once the 6*23 degeneracy is lifted one can easily obtain four allowed sets of (J, 
Am|i) (though they are still degenerate at almost the same point on the sin^ 26*13-823 plane) 
because the relationship between them is given analytically in a completely general setting 

0 . 


B. Resolving power of the (@13, 623) degeneracy by a reactor measurement 

Let us make a semi-quantitative estimate of how powerful the reactor method is for 
resolving the (^13, 6*23) degeneracy.^ ^ For this purpose, we compare in this section the 
difference of the two 6*13 solutions due to the degeneracy with the resolving power of the 


® In the absence of the matter effect, the reason why the closed curve shrinks into a line at the the oscillation 
maximum can be seen as follows: By eliminating S in m, it is easy to show that there are two solutions 
of sin2013 > 0 for given values of P, P and 023 off the oscillation maximum (A 31 ^ tt), whereas there is 
only one solution of sin 20 i 3 > 0 at the oscillation maximum (A31 = tt). Even if we switch on the matter 
effect, one can easily show by using the approximate formula in |l 8 | that the same argument hol(j 
^ The possibility of resolving the ( 0 i 3 , 023) by a reactor experiment was qualitatively mentioned in 
® An alternative way to resolve the ambiguity is to look at Vf, —> Vt channel because the main oscillation 
term in the probability P{i'e z^r) depends up on c?q. Unfortunately, this idea does not appear to be 


explored in detail while it is briefly mentioned in | 2 J, |2 
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reactor experiment. We consider, for simplicity, the special case | A 31 I = tt, i.e., energy tuned 
at the hrst oscillation maximum. The simplest case seems to be indicative of features of 
more generic cases. 

As we saw in the previous section, there are two solutions of 6*13 due to doubling of 6*23 
for a given sin^ 2^23 in each sign of Arug^. Then, we dehne the fractional difference due to 
the degeneracy 

(5de(sin^26'i3) 
sin^ 26*13 


It is to be compared with 5re(sin^ 26 * 13 )/sin^ 26*13 of the reactor experiment, where 
5re(sin^ 26 * 13 ) denotes the experimental uncertainty estimated in Sec. Ill, i.e., 0.043 or 0.018. 
In Fig. a) we plot the normalized error ^^(sin^ 26 * 13 )/sin^ 26*13 which is expected to be 
achieved in the reactor experiment described in Sec. III. We restrict ourselves to the analy¬ 
sis with 1 degree of freedom, because we expect that the JHF phase I will provide us accurate 
information on Arugg by the time when the issue is really focused on the degeneracy in the 
JHF phase II. The fractional difference m can be computed from the relation j 2 ^ 


■ 9 ■ 9 9n / Aruni \ ^ tan^ (aL/2) 

sm^ 20ig = sm^ 2^13 tan^ ^23 + w— ^ , L 

{aL/ny 

X [1 — (aL/vr)^] sin^ 26*12 (l — tan^ 6*23) , 


(13) 


and the result for Jde(sin^ 26 * 13 )/sin^ 26*13 is plotted in Fig. ini(b) as a function of sin^ 6*23 
for two typical values of e. We notice that the fractional differences differ by up to a 
factor of ~ 2 in small sin^ 2023 region between the hrst (023 < vr/4) and the second octant 
(023 > 7 r/ 4 ). For the best £t value of the two mass squared differences Am^i (6.9 x 10“® eV^) 
and |Am|g| (2.5 x 10“^eV^), for which e = Amli/\Aml^\ = 0.028, there is little difference 
between the case with sin^ 20 i 3 = 0.03 and the one with sin^ 20 i 3 = 0.09. In this case 
they are all approximated by the hrst term in m and Jde(sin^ 2013 )/sin^ 2013 depends 
approximately only on 623 , making the analysis easier. On the other hand, if the ratio 
e = Am|g/|Am|g| is much larger than that at the best ht point, then the second term in 
m is not negligible. In Fig. EKb), (5de(sin^ 2013 )/sin^ 20 i 3 is plotted in an extreme case of 
e = 1.9 X lO^'^eV^/l.G x 10“^ eV^ = 0.12, which is allowed at 90% CL (atmospheric) or 
95% CL (solar), with sin^ 20 i 3 = 0.03, 0.06, 0.09. From this, we observe that the suppression 
in the hrst term in m is compensated by the second term for sin^ 20 i 3 = 0.03, i.e., the 
degeneracy is small and therefore resolving the degeneracy is difficult in this case. To clearly 
illustrate the resolving power of the degeneracy by the reactor measurement, assuming the 
best ht value e = 0.028, we plot in Fig. [7| the region where the degeneracy can be lifted 
in the sin^ 20 i 3 -sin^ 2023 plane. It is evident that the reactor measurement will be able to 
resolve the ( 0 i 3 , 623 ) degeneracy in a wide range inside its sensitivity region, in particular 
for 023 in the second octant. 

Quantitative estimation of the signihcance of the fake solution requires detailed analysis 
of accelerator experiments which includes the statistical and systematic errors as well as the 
correlations of errors and the parameter degeneracies, and it will be worked out in future 
communication. 
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VI. MORE ABOUT REACTOR VS. LONG-BASELINE EXPERIMENTS 


The foregoing discussions in the previous section implicitly assume that the sensitivities 
of reactor and LBL experiments with both v and v channels are good enough to detect 
effects of nonzero 6 ^ 13 . However, it need not be true, in particular, in coming 10 years. To 
further illuminate complementary roles played by reactor and LBL experiments, we examine 
their possible mutual relationship including the cases where there is a signal in the former 
but none in the latter experiments, or vice versa. For ease of understanding by the readers, 
we restrict our presentation in this section to a very intuitive level by using a figure. It is, 
of course, possible to make it more precise by deriving inequalities based on the analytic 
approximate formulae H. Throughout this section LBL experiments at the oscillation 
maximum and ^23 = 7 r /4 are assumed. 

If a reactor experiment sees an affirmative evidence for the disappearance in Vf. V^. (the 
case of Reactor Affirmative), it would be possible to determine ^13 up to certain experimental 
errors. In this case, the appearance probability in LBL experiment must fall into the region 
p(^)min ^ Piy) < if the mass hierarchy is known, where the +(—) sign refers to 

Arug^ > 0 < 0 ) and max (min) refers to the maximum (minimum) value of the 

allowed region for P = Piv^ —>■ Ve)-, respectively. (See Fig. |H1) Without the knowledge 
of the mass hierarchy the probability is within the region P(i/)™“ < P(z/) < P(z/)™®'^. 
The similar inequalities are present also for antineutrino appearance channel. In Fig. |Hlwe 
present allowed regions in the cases of > 0 and Arug^ < 0 on a plane spanned by 

PiPfi ^e) and P{u^ —> Pe) by taking two best fit values sin^ 2^13 = 0.08, 0.04 (labeled 
as a, b) as reactor affirmative cases. They are inside the sensitivity region of the reactor 
experiment discussed in section III. We have used the one dimensional analysis (i.e., 
the only parameter is sin^ 2 ^ 13 ) to obtain the allowed regions in Fig. |H1 In doing this we 
have used the same systematic error of 0.8 % and the statistical errors corresponding to 
40 ton-year measurement by the detector considered in section III. For sin^ 2^13 < 0.02, 
the particular reactor experiment would fail (the case of Reactor Negative) but the allowed 
region can be obtained by the same procedure, and presented in Fig. |H1 the region labeled 
as c. We use the same LMA parameters as used earlier for Fig. Inland Fig. HI 

We discuss four cases depending upon the two possibilities of affirmative and negative ev¬ 
idences (denoted as Affirmative, and Negative) in each disappearance and appearance search 
in reactor and long-baseline accelerator experiments, respectively. However, it is convenient 
to organize our discussion by classifying them into two categories, (Reactor Affirmative), 
and (Reactor Negative). 


A. Reactor Affirmative 

We have two alternative cases, the LBL appearance search Affirmative, or Negative. 

LBL Affirmative: 

Implications of affirmative evidence in the appearance search in LBL experiments differ 
depending upon which region the observed appearance probability P{i') falls in: 

( 1 ) < P{v) < or (2) PP^^ < P{v) < P^^^: 

These cases correspond to the two intervals which are given by the projection on the P 
axis of all the shadowed regions (a or b) minus the projection on the P axis of the darker 
shadowed region (a or b) in Fig. |H1 It is remarkable that in these cases not only the sign of 
Am|^ is determined, but also the CP phase 6 is known to be nonvanishing. If F(z/) is in the 
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former region then Am^i is negative and sin <5 is positive, whereas if P{i') is in the latter 
then is positive and sin 5 is negative. 

(3) < P(z/) < 

This case corresponds to the interval which is given by the projection on the P axis of 
the darker shadowed region (a or b) in Fig. |H1 In this case, neither the sign of Am\i nor the 
sign of sin 5 can be determined. 

It may be worth noting that if the reactor determination of di^ is accurate enough, it 
could be advantageous for LBL appearance experiments to run only in the neutrino mode 
(where the cross section is larger than that for antineutrinos by a factor of 2-3) to possibly 
determine the sign of depending upon which region P{i') falls in. 

LBL Negative; 

In principle, it is possible to have no appearance event even though the reactor sees 
evidence for disappearance. This case corresponds to the left edge of the analogous shadowed 
region in the case of sin^ 26*13 — 0.02 in Fig. |Hl i.e., the allowed region with sin^ 20 i 3 ~ 0.02 
for which p™™ on the P axis falls below P = 0.005. In order for this case to occur the 
sensitivity limits P(i^)iimit of the LBL experiment must satisfy pp“ < P(i/)iimit assuming 
our ignorance to the sign of If it occurs that pp“ < P(z/)iimit < P™™, then the sign 

of Am\^ is determined to be minus. 

P(z/)iimit of the JHF experiment in its phase I is estimated to be 3 x 10“^ 3-^ Therefore, 


by using the mixing parameters typical to the LMA solution, the case of LBL Negative cannot 
occur unless the sensitivity of the reactor experiment becomes sin^ 26*13 ~ 0.01. However, in 
the intermediate stage of the JHF experiment, where P(z/)iimit is larger than 3 x 10“^, this 
situation may occur. 


B. Reactor Negative 

If the reactor experiment does not see disappearance of one obtains the bound 
^13 < We have again two alternative cases, the LBL appearance search Affirmative, or 
Negative. 

LBL Affirmative: 

If a LBL experiment measures the oscillation probability P{v)- Then, for a given value 
of P{i') the allowed region of sin 26*13 is given by sin 20 “™ < sin 26*13 < sin 26*™®'^ if the sign 
of Am\i is known, and by sin 20™“ < sin 20 i 3 < sin20™®'^ otherwise. We denote below 
the maximum and the minimum values of 0 i 3 collectively as 0max and 0min, respectively. In 
Fig- in the region bounded by sin 20 “™ and sin 20 ™®'^ (sin 20 “™ and sin 20 ™®'^) are indicated 
as a region bounded by the solid (dashed) black line for a given value of S 23 . 

Then, there are two possibilities which we discuss one by one: 

(i) 6 * 1 ^'" > 0 max: In this case no additional information is obtained by nonobservation of 
disappearance of Pe in reactor experiment. 

(ii) 0min < 013" < 0max: In tliis case we have a nontrivial constraint 0min < 0 i 3 < 0 i 3 ^- 

LBL Negative: 

In this case, we obtain the upper bound on 0 i 3 , which however depends on the assumed 
values of 5 and the sign of Am|^. A J-independent bound can also be derived: 0 i 3 < 


® The sensitivity limit of sin^ 20i3 quoted in sin^ 20i3 < 6 x 10 obtained by using one-mass scale 
approximation (e <C 1) may be translated into this limit for P{v). 
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min[6'^^,6'max]. 

VII. DISCUSSION AND CONCLUSION 

In this paper, we have explored in detail the possibility of measuring sin^ 29i3 using re¬ 
actor neutrinos. We stressed that this measurement is free from the problem of parameter 
degeneracies from which accelerator appearance experiments suffer, and that the reactor 
measurement is complementary to accelerator experiments. We have shown that the sensi¬ 
tivity to sin^ 2013 ~ 0.02 (0.05) is obtained with a 24.3 GWth reactor with identical detectors 
at near and far distances and with data size of 40 (10) ton-year assuming that the relative 
systematic error is 0.8% (2%) for the total number of events. In particular, if the relative 
systematic error is 0.8 %, the error in sin^ 20 i 3 is 0.018 which is smaller than the uncer¬ 
tainty due to the combined (intrinsic and hierarchical) parameter degeneracies expected in 
accelerator experiments. We also have shown that the reactor measurement can resolve the 
degeneracy in 023 ^ 7r/2 — 623 and determine whether 623 is smaller or larger than 7r/4 if 
sin^ 2013 and cos^ 2023 are relatively large. 

We have taken 2 % and 0.8 % as the reference values for the relative systematic error for 
the total number of events. 2 % is exactly the same figure as the Bugey experiment while 
0.8% is what we naively expect in the case we have two identical detectors, near and far, 
which are similar to that of the CHOOZ experiment. It is also technically possible to dig 
a 200 m depth shaft hole with diameter wide enough to place a CHOOZ-like detector in. 
Therefore, the discussions in this paper are realistic. We hope the present paper stimulates 
interest of the community in reactor measurements of 0 i 3 . 
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Bugey 

absolute normalization 

relative normalization 

relative / absolute 

flux 

2.8% 

0.0% 

0 

number of protons 

1.9% 

0.6% 

0.32 

solid angle 

0.5% 

0.5% 

1 

detection efficiency 

3.5% 

1.7% 

0.49 

total 

4.9% 

2.0% 



CHOOZ-like 

absolute normalization 

relative normalization (expected) 

relative / absolute 

flux 

2.1% 

0.0% 

0 

number of protons 

0.8% 

0.3% 

0.38 

detection efficiency 

1.5% 

0.7% 

0.47 

total 

2.7% 

0.8% 


for bins 

8.1% 

2.4% 



TABLE I: Systematic errors in the Bugey and the CHOOZ-like experiments. Relative errors in the 
CHOOZ-like experiment are expectation with the same reduction rates of errors as those of Bugey. 
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CHOOZ, 2d.o.f. - 

Ogys = 2%, 10t»yr, 2d.o.f. . 

Ogys = 2%, 10t*yr, Id.o.f.- 

Ggys = 0.8%, 40t»yr, 2d.o.f. 

Ogys = 0.8%, 40t*yr, Id.o.f. - 
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sin2 29i3 


FIG. 1; Shown are the 90% CL exclusion limits on sin^ 20i3 which can be placed by the reactor 
measurement as described in Sec. III. From the left to right, the dash-dotted and the thin-dotted 
(the long-dashed and short-dashed) lines are based on analyses with 1 and 2 degrees of freedom 
(see the text), respectively for (Tsys=0.8%, 40t-yr {asys=2%, lOt-yr). The solid line is the CHOOZ 
result, and the 90 % CL interval 1.6 x 10“^ eV^ < < 3.9 x 10“^ eV^ of the Super-Kamiokande 

atmospheric neutrino data is shown as a shaded strip. 
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FIG. 2: Shown is the accuracy of determination of sin^ 2^13 at 90 % CL for the case of positive 
evidence based on analysis with 1 degree of freedom, = 2.7. Figures (a) and (b) are for 

crsys=2%, 10 t-yr, and cJsys=0.8%, 40 t-yr, respectively. The lines correspond to the best fit values 
of sin^ 2013, from left to right, 0.05 to 0.08 in the unit of 0.01 in Fig. lU^a), and 0.02 to 0.08 in the 
unit of 0.01 in Fig. (Jl^b). The reference value of |Am 3 i*’'^^*^| is taken to be 2.5 x 10“^ eV^, which is 
indicated by a gray line. 
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FIG. 3: Depicted in the sin^ 2613-823 plane are the contours determined by arbitrarily given 
values of the appearance probabilities P = ^ I'e) = O-Ol P = P{v^ i^e) = 0.015 

with E/L off the oscillation maximum (IA 31 I 7 ^ vr) at the JHF experiment. Here, S 23 = sin^023- 
The solid and the dashed lines correspond to positive and negative respectively. The dash- 

dotted lines represent the boundary of the region 0.36 < S 23 — 0-64 which is presently allowed 
by the atmospheric neutrino data, 0.92 < sin^ 2,623 < 1- As indicated in the figure, there are 
four solutions for each S 23 ) altogether there are eight solutions as denoted by blobs for any 
values of 623 7 ^ vr/d. The oscillation parameters are taken as follows: Am|^ = 2.5 x 10“^eV^, 
Am|^ = 6.9 X 10“^eV^, tan^ 6*12 = 0.38. The Earth density is taken to be p=2.3 g/cm^. 
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FIG. 4: The allowed regions are shown in the sin^ 20 i 3 -s |3 plane determined with a given value of 
P = —> t'e) alone (in this case P = 0.025), or P = P{i>^ —> Ve) alone (in this case P = 0.035) 

at the oscillation maximum IA 31 I = vr of the JHF experiment. Each allowed region is the area 
bounded by the black solid (for > 0 with P only), the black dashed (for Am|^ < 0 with P 

only), the gray solid (for Am|^ > 0 with P only), the gray dashed (for Am|^ < 0 with P only), 
respectively, where the line with a definite value of the CP phase 6 sweeps out each region as 6 
varies from 0 to 27r. The oscillation parameters and the Earth density are the same as those in 

Fig.H 
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P=0.025, P=0.035, Am§i >0- 

P=0.025, P=0.035, Am§i <0- 



FIG. 5; The allowed region in the sin^ 20 i 3 -s |3 plane becomes a line when both P{Vf^ fe) and 
^e) are given (in this case P{Vfjt Ve) = 0.025, P{i>f_i i'e) = 0.035) at the oscillation 
maximum (IA 31 I = tt, E = 0.6 GeV for the JHF experiment), as indicated in the figure. The solid 
and the dashed lines are for > 0 and Am^i < 0 cases, respectively. Assuming 623 / vr/d, 

two solutions of (sin^ 20 i 3 , S 23 ) are plotted; In this figure sin^ 26*23 is taken as 0.92. It is assumed 
arbitrarily that the solution of 6*23 in the first octant (023 < vr/d) is the genuine one, while the one 
in the second octant (023 > vr/d) with primes is the fake one. Superimposed in the figure as a 
shaded region is the anticipated error in the reactor measurement of 0 i 3 estimated in Sec. III. If 
the error (5re(sin^ 20 i 3 ) is smaller than the difference (Ide(sin^ 20 i 3 ) = | sin^ 20^3 — sin^ 20 i 3 | due to 
the degeneracy, then the reactor experiment may be able to resolve it. 


22 
















eggUjS / (S^022U!S)®^§ 



FIG. 6 : (a) The normalized error at 90 % CL in the reactor measurement of @13 is given 

for (Tsys=2%, lOt-yr (d.o.f.=l, (Jre(sin^ 26 * 13 ) = 0.043) and for (Tsys=0.8%, 40t-yr (d.o.f.=l, 
(5re(sin^ 2 ^ 13 ) = 0.018), respectively. Notice that the degrees of freedom becomes 1 once the value 
of is known from JHF. 

(b) The fractional difference (Jde(sin^ 26 * 13 )/ sin^ 2^13 due to the degeneracy is plotted as a function of 
sin^ 26 * 23 . Here, (5de(sin^ 2 ^ 13 ) = | sin^ 20 ^ 3 —sin^ 20 i 3 | stands for the difference between the true solu¬ 
tion sin^ 2013 and the fake one sin^ 20 ^ 3 , and e = Amli/IAmlj^l; e = 6.9x 10“® eV^/2.5x 10“^ eV^ = 
0.028 is for the best fit and an extreme case with e = 1.9 x 10“^ eV^/1.6 x 10“^ eV^ = 0.12, which is 
allowed at 90% CL (atmospheric) or 95% CL (solar), is also shown for illustration. The horizontal 
axis is suitably defined so that it is linear in sin^ 2023 > where the left half is for 623 < 7 r /4 whereas 
the right half is for 623 > 7 r/ 4 . The solar mixing angle is taken as tan^ 0 i 2 = 0.38. sin2 2023 > 0.92 
has to be satisfied due to the constraint from the Super-Kamiokande atmospheric neutrino data. 
If the value of cos^ 2023 is large enough, the value of (5de(sin^ 2013 )/ sin^ 20 i 3 increases and lies out¬ 
side of the normalized error of the reactor experiment, then the reactor result may resolve the 023 
ambiguity. 
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FIG. 7: The shadowed area stands for the region in which (5re(sin^ 2 ^ 13 ) < 5de(sin^ 2 ^ 13 ) is satisfied 
for (Tsys=0.8%, 40t-yr, d.o.f.=l and for the best fit values of the solar and atmospheric oscillation 
parameters. In this shadowed region, the (0i3 > ^ 23 ) degeneracy may be solved. The vertical axis is 
the same as the horizontal axis of Fig. Hb). 
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FIG. 8 : Predicted allowed regions are depicted in the P-P plane for the JHF experiment at the 
oscillation maximum after an affirmative (a negative) result of the reactor experiment is obtained, 
where P = Pi^^ r'e) and P = P{i>ij, r'e) are the appearance probabilities, and 623 = vr/d 
is assumed. The cases a, b, c correspond to sin^ 26*13 = 0.08 ± 0.018, sin^ 2^13 = 0.04 ± 0.018, 
sin^ 2013 < 0.019, respectively. The regions bounded by the solid lines and the dashed lines are 
for the normal hierarchy {Am^i > 0) and the inverted hierarchy (Am|^ < 0), respectively. Each 
region predicts the maximum and the minimum (P™“) values of P for each hierarchy (+ 

for the normal and — for the inverted hierarchy), although of the region c are zero. 
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